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ABSTRACT 

Atmospheric  turbulence  has  been  shown  to  have 
measurable  effects  on  the  angle-of-arrival  (AO A)  of  elec¬ 
tromagnetic  wave  beams,  but  these  effects  are  on  the 
order  of  a  few  microradians,  which  is  a  negligible  level 
for  most  applications.  In  this  paper,  we  present  a  theory 
describing  this  phenomenon  and  compare  this  theory  to 
one-way  measurements  made  over  a  3.5  km  path  and  two- 
way  measurements  made  over  a  25  km  path  using  an  X- 
band  interferometric  radar. 

1.  INTRODUCTION 

The  effects  of  atmospheric  turbulence  on  radar  angle 
errors  are  generally  small,  but  in  those  cases  where  a 
radar  is  required  to  guide  a  missile  without  a  seeker 
(command  guide),  a  few  microradians  of  error  may  be 
enough  to  cause  the  missile  to  miss  its  target.  In  this 
paper  we  compare  two  methods  of  calculating  these  an¬ 
gle-of-arrival  (AO A)  errors  and  show  that  the  geometrical 
optics  approach  derived  by  Chumside  and  Lataitis  [1],  the 
first  method,  gives  larger  angular  errors  than  those  deter¬ 
mined  by  experiment.  This  method  is  characterized  by 
perfectly  collimated  and  focused  beams  and  ideal  plane 
waves.  The  second  method  involves  the  adaptation  of  the 
approach  described  in  [1]  to  a  physical  optics  model  in 
which  the  beams  have  a  Gaussian  profile.  This  method  of 
calculation  gives  results  reasonably  close  to  experimental 
values.  We  determine  the  AO  A  using  both  approaches 
and  compare  the  results  obtained  to  an  experiment  con¬ 
ducted  using  an  interferometric  radar  capable  of  measur¬ 
ing  the  small  angles  expected  for  this  scenario. 

2.  THEORY 

In  Reference  [I],  the  authors  have  derived  the  AO  A  of 
a  focused  beam  in  the  geometrical  optics  limit  for  a  one¬ 
way  path  and  have  shown  that  their  result  reduces  to  those 
obtained  for  plane  and  spherical  waves  in  the  limits  of 
infinite  and  zero  focal  lengths,  respectively.  They  also 
show  that  the  AOA  approaches  infinity  in  a  predictable 
way  as  range  approaches  the  focal  length  of  the  focusing 
optic.  They  used  the  simple  model  shown  in  Figure  1  as 


the  basis  of  their  calculations.  The  tilt  angle  da  is  given 
by 

da  =  /S.n{z)dz !  w{L),  (1) 

so  that  the  total  tilt  angle  over  path  L  is 

I  ^ 

«=— —  J  Kn{z)dz.  .2') 

w{L)  0  ^  ^ 

Using  this  result,  and  assuming  that  there  is  no  average 
gradient  of  refractive  index,  Chumside  and  Lataitis  derive 
an  expression  for  the  one-way  AOA  based  on  geometrical 
optics. 


Figure  1.  A  thin  atmospheric  layer  with  varying  index 
causes  beam  steering 

Extending  the  arguments  leading  to  Equation  (1)  to 
the  two-way  case,  we  find  that  the  tilt  back  at  the  trans¬ 
mitter  is  given  by 


da  - 


AUf  (z)  -h  AHj.  (z) 

w,(0) 


(3) 


where  Ant  and  An^  are  the  refractive  index  differences 
across  the  transmitted  and  reflected  beams,  respectively, 
and  Wi-(O)  is  the  diameter  of  the  reflected  beam  at  the 
transmitter.  Summing  contributions  along  the  entire  path 
as  for  Equation  (2)  gives  a  total  tilt  angle  of 
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Just  as  Equation  (2)  is  the  basis  for  derivation  of  the  ex¬ 
pressions  for  AOA  for  one-way  propagation,  Equation  (4) 
is  used  to  derive  expressions  for  two-way  propagation. 
The  details  of  these  derivations  are  given  in  [1]  and  will 
not  be  repeated  here. 

In  deriving  the  AOA  for  the  physical  optics  case,  we 
proceed  exactly  as  for  the  geometrical  optics  case  using  the 
procedure  given  in  [1].  For  one-way  transmission,  we  have 
shown  [2]  that  the  variance  of  the  AOA  measured  at  the 
target  is 


)V'’dz  ,5) 

w  (L)  ’  ^  ^ 

where  w(L)  and  w(zi)  are  the  beamwidths  measured  at 
ranges  L  and  z,  respectively  and  Cn^  is  the  index  of  refrac¬ 
tion  structure  parameter.  In  the  physical  optics  (PO)  case, 
the  beam  has  a  Gaussian  profile  and  its  width  varies  as 
[3,4]: 


Substituting  Equation  (6)  into  Equation  (5),  we  get  for  the 
AOA  variance 
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dz. 


(8) 


where  all  parameters  have  been  defined  previously.  In 
the  derivation  of  Equation  (8)  we  assume  that  all  of  the 
radiation  from  the  transmitter  is  collected. 


For  the  two-way  case,  we  simply  substitute  Equation 
(6)  for  the  Gaussian  beam  profile  into  the  two-way  AOA 
expression  (7).  The  resulting  equation  must  be  solved 
numerically  as  must  Equation  (8). 

To  obtain  the  geometrical  optics  result  we  use  an 
equation  derived  in  [1]  that  is  based  on  the  geometrical 
optics  scenario  of  a  collimated  transmitter  beam  illumi¬ 
nating  a  reflector  with  a  negative  focal  length.  We  choose 
the  focal  length  to  be  negative  because  most  targets  of 
interest  will  be  convex.  This  result  is  Equation  (31)  of 
Reference  [1]  and  gives  the  variance  of  the  AOA,  as: 


w(z)  =  2wq 


(6) 


where  X  is  wavelength  and  wq  is  the  1/e  beamwidth  at  the 
transmitter. 

For  a  reflected  beam,  this  approach  is  extended  to  2- 
way  transmission  using  the  method  described  in  [1].  The 
result  for  the  AOA  variance  back  at  the  radar  receiver  is: 


=  imcIld;^‘\  (9) 

where  L  is  slant  range,  Dt  is  antenna  diameter,  and  the 
parameter  Cn^  is  the  refractive  index  structure  parameter. 
In  deriving  this  result,  we  have  used  the  approximation 
L»fi.,  where  f^  is  the  (negative)  focal  length  of  the  reflec¬ 
tor.  We  do  not  have  to  know  the  value  of  the  focal  length 
because  if  we  use  this  approximation,  the  entire  bracketed 
term  in  Equation  (31)  of  Reference  [1]  reduces  to  2.40. 


[h(z)r^^+k(z)r'^ 

1 

^k(z)  +  W,(z)f'^  ^ 
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where  Wi-(O)  is  !/2(l/e)  times  the  diameter  of  the  radar 
antenna,  Wt(z)  is  the  diameter  of  the  transmitter  beam,  and 
Wi-(z)  is  the  diameter  of  the  reflected  beam.  These  latter 
parameters  are  determined  by  substituting  E2(l/e)  times 
the  diameter  of  the  transmitter  and  the  reflector,  respec¬ 
tively.  The  reflector  is  assumed  to  be  a  circular  mirror 
normal  to  the  direction  of  propagation  of  the  transmitter 
beam  for  these  calculations. 

To  determine  the  one-way  AOA  using  the  physical 
optics  approach,  we  must  integrate  Equation  (5)  numeri¬ 
cally  using  the  beam  diameter  given  by  Equation  (6). 


3.  EXPERIMENT 

Measurements  of  AOA  were  made  on  two  separate 
occasions  at  a  location  in  Brea,  CA.  The  first  series  of 
measurements  was  made  over  a  path  length  of  25  km 
using  a  passive  reflector  located  on  a  hillside  such  that  the 
elevation  angle  was  about  four  degrees.  Results  from  this 
series  of  measurements  are  presented  in  this  paper.  The 
second  series  of  measurements  was  made  in  the  same 
locality  over  a  range  of  3.5  km  but  with  active  repeaters 
instead  of  reflectors.  Since  these  latter  measurements 
were  made  with  repeaters,  they  are  considered  to  repre¬ 
sent  a  one-way  path,  so  that  Equation  (8)  is  applicable. 
These  latter  measurements  also  gave  very  useful  values  of 
Cn^,  which  were  used  in  the  calculations  presented  herein, 
both  for  the  one-way  and  the  two-way  paths. 

Figure  2  shows  the  experimental  arrangement  used 
for  the  long-range  measurements  and  Figure  3  shows  the 
setup  for  the  short-range  measurements  and  those  used  to 


determine  Cn^.  Both  of  these  experiments  use  interfer¬ 
ometric  radars  for  measurements  of  the  very  small  AO  As 
characteristic  of  microwave  frequencies.  For  determina¬ 
tion  of  Cn^,  the  fluctuations  from  a  single  channel  were 
used  and  Cn^  was  calculated  from  the  log  amplitude  vari¬ 
ance  using  the  relation 

cr^  (10) 

where  k  is  wavenumber  InlX. 


Figure  2.  Diagram  of  the  interferometric  radar  used  for 
long-range  AO  A  measurements. 


Figure  3.  Diagram  of  the  interferometric  radar  used  for 
short-range  AO  A  measurements  and  determination  of  Cn^. 
These  measurements  were  made  at  the  Technovative 
Applications  facility  in  Brea,  CA. 

Figure  4  is  a  photograph  of  the  repeaters  used  for 
short-range  one-way  measurements.  The  antenna  aper¬ 
tures  are  12.3  X  9.1  cm.  Figure  5  is  a  photograph  of  the 
interferometric  radar.  The  transmitter  is  in  the  center  and 
the  receivers  are  at  the  ends  of  the  boom  with  a  spacing  of 
5  m.  The  transmitter  frequency  was  9523  MHz  and  the 
ntenna  diameter  was  56  cm. 


Figure  4.  Photograph  of  the  repeaters  used  for  short- 
range  AO  A  and  Cn^  measurements. 


Figure  5.  Photograph  of  the  interferometric  radar  used 
for  long-range  AO  A  measurements. 

4.  RESULTS 

It  was  first  necessary  to  determine  the  values  of  Cn^ 
for  the  short-range  experiment.  Calculations  of  AOA 
were  then  made  using  Equation  (8)  with  these  values  and 
compared  to  measured  results.  Other  measurements  of 
Cn^  made  at  millimeter  waves  [5]  have  shown  that  this 
parameter  can  be  as  high  as  5.5  X  10  m'^^^  under  hot, 

humid  conditions.  These  parameters  are  much  smaller  in 
the  visible  and  infrared  range  and  partially  explain  the 
result  that  turbulent  fluctuations  are  still  observable  in  the 
microwave  bands  even  though  theory  predicts  that  they 
should  decrease  as  where  f  is  frequency. 


We  have  solved  Equation  (8)  numerically  for  ranges 
between  500  and  5000  meters,  a  transmitter  aperture  of 
0.1  m,  a  frequency  of  9.5  GHZ,  and  values  of  Cn^  varying 
from  1.4  X  10'^^  to  6.8  X  10'^^  and  the  results  are 
shown  in  Figure  6. 


used  1.0  X  10'^^  m'^^^  for  this  parameter,  which  is  more 
consistent  with  reality  based  on  Figure  6,  we  would  have 
obtained  a  standard  deviation  of  9  microradians,  in  good 
agreement  with  the  results  shown  in  Figure  9. 


Figure  6.  Calculated  AOA  variance  as  a  function  of  range 
for  the  one-way  path  for  the  conditions  described  in  the 
text. 


Table  I  summarizes  these  results  for  the  3.5  km  range 
used  in  the  measurements.  These  are  the  physical  optics 
(PO)  results  shown  in  the  table.  Calculation  of  the  geo¬ 
metrical  optics  (GO)  numbers  were  made  using  the  equa¬ 
tions  of  Reference  [1]  and  the  values  of  Cn^  were  deter¬ 
mined  using  Equation  (1).  This  level  of  agreement  is 
considered  good  for  this  type  of  experiment  for  both  sets 
of  calculations,  but  the  results  show  that  the  PO  formal¬ 
ism  gives  better  agreement  with  experiment. 

We  have  used  Equation  (7)  to  calculate  the  AOA  ex¬ 
pected  for  the  25  km  two-way  path.  Although  simultane¬ 
ous  measurements  of  Cn^  were  not  made,  the  value  of  this 
parameter  corresponding  to  Figure  5  above  was  consid¬ 
ered  typical  and  was  used  in  these  calculations.  Figure  7 
shows  the  results  obtained  for  values  of  the  diameter  of 
the  remote  reflector  of  0.2,  0.4,  0.6,  and  0.8  meters.  For 
0.6  meters  at  25  km  range,  these  curves  show  that  the 
AOA  variance  is  7.3  X  10'^^  rad^  corresponding  to  a  stan¬ 
dard  deviation  of  27  microradians.  The  value  of  Cn^  used 
in  these  calculations  was  3.3  X  10'^^  and  the  trans¬ 
mitter  antenna  diameter  used  was  0.6  m. 

Figure  9  is  a  scatter  plot  of  cross  range  errors  meas¬ 
ured  with  the  interferometric  radar.  This  figure  shows 
that  the  standard  deviation  of  these  errors  is  about  10 
microradians.  This  agreement  between  measured  and 
calculated  values  of  AOA  is  considered  to  be  fair.  In  the 
next  section,  we  will  discuss  the  possible  reasons  for  this 
level  of  disagreement.  Note  that  the  parameter  was 
not  measured  for  the  long-range  results.  We  used  the 
value  3.3  X  10'^^  for  these  calculations.  If  we  had 


RANGE  IN  KILOMETERS 

Figure  7.  Results  of  calculating  AO  for  two-way  paths 
varying  from  5  to  50  km  for  various  receiver  reflector 
diameters.  The  parameter  of  interest  for  this  experiment 
are  25  km  range  and  0.6  m  reflector  diameter. 


CROSS  RANGE  AZIMUTH  ERROR  [cm] 


Figure  8.  Scatter  plot  of  AOA  measurements  made  with  a 
0.6m  transmitter  and  reflector  over  a  25  km  path.  The 
measured  standard  deviation  is  10  microradians 

It  is  instructive  to  compare  the  results  depicted  above 
to  those  obtained  by  using  the  one-way  equation 
(jf  =  2.92LCId~^'^  ,  where  Dj-  is  receiver  diameter,  and 
the  two-way  equation  as  derived  in  [1].  Using  the  same 
parameters,  we  get  an  AOA  standard  deviation  of  63 
microradians  for  the  one-way  case  and  262  microradians 
for  the  two-way  case.  Even  though  the  analysis  given 


above  gives  values  of  AOA  that  are  consistently  higher 
than  measured  values,  they  represent  a  considerable  im¬ 


provement  over  the  results  obtained  by  using  the  earlier 
accepted  techniques  for  calculating  AOA. 


Table  L  Results  Of  AOA  Measurements  Over  A  One-Way,  3.5  Km  Path. 


AOA  Standard  Deviation  in  Microradians 

CTinnF^^ 

5.8X10-'^ 

2.0X10-'^ 

1.4X10''^ 

2.6X10-*^ 

6.3X10-'^ 

GO  Calculated 

5.9 

3.5 

2.9 

3.9 

6.1 

PO  Calculated 

1.6 

0.9 

0.8 

1.1 

1.7 

Measured 

0.9 

0.5 

0.3 

1.2 

2.8 

5.  CONCLUSIONS 

We  have  developed  a  method  of  calculating  radar 
AOA  based  on  the  physical  optics  principle  of  Gaussian 
beam  profiles.  In  this  approach,  the  radar  beams  are  al¬ 
lowed  to  diverge  normally  instead  of  being  perfectly 
collimated  as  in  the  geometrical  optics  case. 

We  find  that  our  theory  predicts  AOAs  that  are  con¬ 
sistently  slightly  higher  than  measured  values.  Part  of 
this  error  can  be  ascribed  to  the  slight  elevation  of  the 
propagation  path  and  the  resulting  slight  reduction  in  Cn^ 
as  a  function  of  altitude  for  which  our  theory  does  not 
account.  Other  possible  sources  of  error  are  noise  gener¬ 
ated  in  the  repeater  receivers  and  transmitters  and  noise  in 
the  radar  receiver.  Other  problems  include  repeater  or 
radar  motion  caused  by  wind  or  by  ground  vibrations, 
although  all  of  these  potential  errors  would  tend  to  in¬ 
crease  AOA  standard  deviation  instead  of  reducing  it. 

Despite  this  lack  of  agreement  with  measured  AOA, 
our  results  are  still  better  than  those  obtained  by  using  a 
geometrical  optics  approach.  The  poor  agreement  using 
geometrical  optics  results  from  the  lack  of  diffractive 
beam  spreading  inherent  in  the  assumptions  used  in  deriv¬ 
ing  this  theory.  We  expect  to  continue  measurements  of 
AOA  using  the  apparatus  described  herein  and  will  refine 
our  theory  to  attempt  to  get  better  agreement  between 
theory  and  measurement. 
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“Secure  the  High  Ground” 


Purpose 


•  Precise  radar  applications  require  high  angular 
accuracy 

•  It  is  necessary  to  determine  the  atmospheric 
contribution  to  angular  error 

•  Knowledge  of  the  index  of  refraction  structure 
parameter  0^^  is  necessary  to  quantify  angular 
errors 

•  Measurements  of  angle-of-arrival  are  necessary 
to  support  our  theoretical  basis 
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Procedure 

•  Measure  fluctuations  in  X-band  signal  over 
a  path  of  suitable  length  L 

•  Calculate  using  the  equation 

-  Where  log  amplitude  variance  and  k  is 
wavenumber 

•  Calculate  AOA  using  our  theory 

•  Compare  theory  and  experiment 
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A  Thin  Atmospheric  Layer  With  Varying  Index 

Causes  Beam  Steering 


da  =  hn{z)dz  I  w{L) 

L 


In  deriving  these 
equations,  we  assume  that 
there  is  no  average  gradient  of 
refractive  index 
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dz 


One-Way  Path 


The  variance  of  AOA  fluctuations  for  a  one-way  path  is  then 
given  by: 


w 


(An(z^)An(z2))dz^dz2, 


The  index  gradient  can  be  written  in  terms  of  two  spatial 
coordinates  as 


-  n 


z- 


5 
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Beam  Geometry 
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One-Way  Path  (Continued) 

The  quantity  in  the  angle  brackets  can  be  expanded  and 
rearranged  to  get 

/ {«[Zj  ,1  /  2w(z, )]  -  «[Z2,1/2w(Z2)]}^ 

(A„(.,  M-'z )) = 1 , 2/ + )]  -  »[z„i  /  2M.Z  )1 

\  -  {n[zj  ,1  /  2w(zj )]  -  «[z2,-1/2w(z2)]}^ 

{n[zj,-l/2w(zj)]  -  «[z2,-1/2w(z2)]}^ 

Using  the  equation  for  the  calculation  of  the  distance  between 
two  points  and  and  the  definition  of  the  index  of  refraction 
structure  parameter  Cl={n{r^)-n{r2))^ l{r^ we  find  that  for 
homogeneous,  isotropic  turbulence,  the  above  quantity  becomes 
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One-Way  Path  (Continued) 


(An(z,  )A«(z2  ))  =  Cl  ({(Z2  -  Zj  )^  + 1  / 4[w(z2 )  +  w(z^ 

-[(^2  +l/4[>v(Z2)->v(Zi)f 

Now  assume  that  the  magnitude  of  the  focal  length  of  the  optical 
system  that  determines  these  beamwidths  is  much  larger  than 
the  beam  diameter  so  that  the  quantity  is  a  slowly  varying  function 
of  position  along  the  path.  For  a  radar  beam,  this  condition  is 
easily  met  because  the  focus  is  at  infinity.  We  may  therefore 
consider  the  beam  diameter  to  be  constant  unless  the  separation 
becomes  much  greater  than  the  diameter.  In  this  case,  the  above 
equation  becomes  very  small  because  the  turbulent  fluctuations 
become  decorrelated.  Under  these  assumptions,  w(zi)  «  w(z2) 
and  the  above  integral  may  be  written - - 


One-Way  Path  (Continued) 


Most  of  the  contribution  to  the  integral  comes  from  values  of  Z2 
that  are  very  nearly  equal  to  z^,  so  there  is  little  error  involved  in 
extending  this  integral  from  negative  infinity  to  infinity.  Note  that 
the  above  integral  is  then  of  the  form 


lo  +  -i/r]dydz, 


Integrating  over  Z2  numerically  gives  simply 

2.92CI 


Jo  Jo  = 


w, 


0 
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For  One-Way  Path 

^  )Y^^dz 


a;  =2.92  ^ 


w\L) 

From  the  theory  of  Gaussian  beam  propagation, 


w(z)  =  2wq 

Substituting 


( 

1  +  \^Z  /  ^TWq  jt 


2.92C 


n 


(2wJ'fl +  (/!£/ 


tJ, 


TIW, 


0 


L 

0 


1  + 


\^Z  /  TTWq  j 


dz 


J.  H.  Churnside  and  R.  J.  Lataitis,  “Angle  of  Arrival  Fluctuations  of  a  Reflected  Beam  in 
Atmospheric  Turbulence”,  J.  Opt  Soc.  Am.  4,  July  1987,  1264. 
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For  Two-Way  Path 


da  = 


An^  (z)  +  AUj.  (z) 


w,(0) 


L 

|K(z)  +  A«,(z)}/z 
0 
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For  Two-Way  Path 


Subject  to  several  assumptions,  the  expression  for  a  reflected 
beam  (two-way  path)  becomes: 


V 


5/3  ,  r  (  n15/3 

+  K  (^)  J 

\wt{z)  +  Wr{z)]‘^ 

^-|[w^(z)-w^(z)] 


5/3 


>V 


J 


>dz 


If  we  now  substitute  the  expression  for  Gaussian  beamwidth 
given  above  for  w/z)  and  w/z),  we  obtain  an  equation  that 
Can  be  solved  numerically. 

-  “Secure  the  High  Ground”  - 


Geometrical  Optics  AOA  Theory 


=2.92C^Lw~^^^ 

For  a  one-way  path  and  aperture  w 
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mm 


Transmitter  (center)  and  Interferometric^ 

Receiver 
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One-Way  Measurements 


We  used  repeaters  (one  at  a  time) 


Repeater  Diagram 


3.5  km 
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Determination  of 


•  No  instrumentation  to  measure  directly 

•  Used  signal  from  one  of  the  receivers  and 
the  equation 

2  A  "2  1  2  /  7/6t'11/6 

Gy  =  0.31C„  A:  L  , 

•Values  measured  were  somewhat  lower  than 
those  generally  observed  for  microwave 
systems,  which  is  inconsistent  with  earlier  work* 

*R.  J.  Hill,  R.  A.  Bohlander,  S.  F.  Clifford,  R.  W.  McMillan,  J.  T.  Priestley, 

W.  P.  Schoenfeld,  “Turbulence-Induced  Millimeter-Wave  Scintillation  Compared 
with  Micrometeorological  Measurements”,  IEEE  Trans.  Geosciences  and  Remote 
Sensing,  26,  May  1988,  330. 
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3:06  PM 


9:4SPM 


6:19  PM 


1 1 :53  AM  The  next  day 


Calculated  One-Way  AOA  Standard  Deviation 


RANGE  IN  KILOMETERS 
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Histograms  of  AOA 
Colors  Correspond  to  Previous  Slide 


Broadside 


30^  Off  Broadside 


Tabular  Results 

Geometrical  Optics  (GO),  Physical  Optics  (PO) 


AOA  Standard  Deviation  in  Microradians 

in  m'2/3 

5.8X10- 

13 

2.0X10-13 

1.4X10-13 

2.6X10-13 

6.3X10-13 

GO 

Calculated 

5.9 

3.5 

2.9 

3.9 

6.1 

PO 

Calculated 

1.6 

0.9 

0.8 

1.1 

1.7 

Measured 

0.9 

0.5 

0.3 

1.2 

2.8 

“Secure  the  High  Ground” 


mmi 


Scatter  Plot  of  AOA  Measurements 


25  km  Range 


CROSS  RANGE  AZIMUTH  ERROR  [cm] 
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Calculated  2-Way  AOA  Versus  Range 
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Conclusions 

•  We  have  measured  using  x-band 
signals  over  3.5  km  path 

•  Correlation  with  AOA  gives  good  results 

•  We  used  results  to  calculate  2-way 
AOA  over  25  km  path 

-  Good  agreement  with  measured  results 

-  Geometrical  optics  calculations  also  give  good 
results 
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